Bright-like soliton solution in quasi-one-dimensional BEC in third order on interaction radius 
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Nonlinear Schrodinger equations and corresponding quantum hydrodynamic (QHD) equations are 
widely used in studying ultracold boson-fermion mixtures and superconductors. In this article, we show 
that a more exact account of interaction in Bose-Einstein condensate (BEC), in comparison with the Gross- 
Pitaevskii (GP) approximation, leads to the existence of a new type of solitons. We use a set of QHD 
equations in the third order by the interaction radius (TOIR), which corresponds to the GP equation in a 
first order by the interaction radius. The solution for the soliton in a form of expression for the particle 
concentration is obtained analytically. The conditions of existence of the soliton are studied. It is shown 
what solution exists if the interaction between the particles is repulsive. Particle concentration of order 
of 10^2-10" cm has been achieved experimentally for the BEC, the solution exists if the scattering 
length is of the order of 1 ^va, which can be reached using the Feshbach resonance. It is one of the limit 
case of existence of new solution. The corresponding scattering length decrease with the increasing of 
concentration of particles. The investigation of effects in the TOIR approximation gives a more detail 
information on interaction potentials between the atoms and can be used for a more detail investigation 
into the potential structure. 



I. INTRODUCTION 



Nonlinear structures, solitons, and vortices have been ac- 
tively studied in atomic Bose-Einstein condensate (BEC) 
and boson-fermion mixtures situated in magnetic and op- 
tical traps lH]]- isj]- Solitons exist in inquest of nonlin- 
earity of interactionally conditioned terms. These nonlin- 
ear terms compensate the dispersion emerging particularly 
as a consequence of a free motion of quantum particles. 
Solitons in the BEC are described via the Gross-Pitaevskii 
(GP) equation |6] which has form of one-particle non- 
linear Schrodinger equation (NLSE). The NLSE plays an 
important role when describing the dynamics of various 
physical systems; name just a few degenerated chargeless 
bosons and fermions as well as superconductors. A more 
detailed account of the interaction in comparison with the 
GP equation leads to the appearance of additional terms 
both in the GP equation and in the corresponding quantum 
hydrodynamics (QHD) equations. A similar generalization 
is acquired for ultracold boson-fermion mixtures [7J. Such 
generalization leads to that the term depending on the spa- 
tial derivatives of the concentration (modulus quadrate of 
wave -function in the medium) appears in the NLSE. Differ- 
ent authors iH], |@] suggested the NLSEs for the descrip- 
tion of the BEC dynamic taking into account the terms de- 
pending on high degrees of concentration. The occurrence 
of additional terms in the GP equation leads to varyin g th e 
characteristics of wave perturbations [10], solitons [lllll . 
vortices, dispersion shock waves, and can also lead to new 
types of solutions. The last point is considered in this pa- 
per. 



The possibility to derive the GP equation from a micro- 
scopic many-particle Schrodinger equation (MPSE) is sub- 
stantiated in Iil2i1 . The method of direct derivation of the 
GP equation from the MPSE was suggested in Ref. 10] ■ It 
was made by means of QHD method Ijj, H. It is well- 
known that the GP equation can be present in the form of 
hydrodynamic equations [6] 



atn(r,t) + a"(n(r,t)t;"(r,t)) = 



(1) 



and 



mn(r, t)dtV°'{r, t) + -mn(r, t)d°'v'^{r, t) 
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+ gn{r, t)a"n(r, t) = -n(r, t)5"V;.t(r, t), (2) 



where 
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drU{r) 



and n(r, t) is the concentration of particles and f "(r, t) is 
the velocity field. The quantity A is the Laplace operator. 

For dilute gases the quantity g can be express via scat- 
tering length by formula 

9 = , 

m 

where a is the scattering length. 

In this paper we use set of equations derived with the 
QHD method. There are different methods for obtaining 
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equations describing the BEC evolution. For example, in 
Ref. 1 15], equation for BEC evolution was derived in the 
framework of nonequilibrium Thermo Field Dynamics. In 
the QHD method the system of equation is appeared di- 
rectly from many-particle Schrodinger equation. The first 
step of derivation is the definition of concentration of par- 
ticles in three dimensional physical space. Differentiation 
of concentration with respect to time and applying of the 
Schrodinger equation leads to continuity equation, a cur- 
rent of density is arisen there. Next step of derivation is 
differentiation of the current density. In this way we obtain 
a momentum balance equation, in another terms the Euler 
equation. A force field exists in the obtained Euler equa- 
tion. For neutral particles with the short-range interaction 
the force field could be present in the form of a expansion 
in a series. In this case, the GP equation emerges when we 
take into account the first member of decomposition of the 
force field by the interaction radius. The next nonzero term 
appears in the third order by the interaction radius (TOIR). 

Different types of solitons occur in the BEC and boson- 
fermion mixtures. If the interaction between the bosons is 
repulsive a > 0, dark solitons that are the regions with 
a lowered concentration of the particles can propagate in 
the BEC ill], 113], iQ. Bright solitons, i.e., solitons 
of compression, can exist in the system of Bose particles 
coupled by attractive forces in quasi-one-dimensional (ID) 
traps lfl9n . lEoll . Gap solitons manifest themselves in 
periodic structures, particularly, the existens of gap bright 
solitons are found experimentally in the system of bosons 
with a > ll2lll . Solitons of compression occur in 
boson-fermion mixtures if repulsive forces a^f, > act 
between the bosons, while the interaction between bosons 
and fermions is attractive with force a^^ < [22]. In the 
work [11] authors obtain a change of the form of the well- 
known bright soliton due to TOIR terms. The bright soliton 
solution arise from GP equation. More detailed account of 
interaction with accuracy to TOIR leads to change of form 
of bright soliton. 

In this paper, we report about the existence of a new soli- 
ton solution in a one dimensional (ID) BEC. This solution 
appears when we account the interaction accurate to the 
TOIR. To obtain this soliton we solve the set of the QHD 
equations by the perturbative method suggested Washimi 
et al. which is widely used in the plasma physics, 

see for example [|24ll . [25]. The obtained solution is the 
soliton of compression; it exists under the condition a > 0, 
i.e., in the case of repulsion between the particles. The ex- 
istence of such a solution can be conditioned by higher spa- 
tial concentration derivatives in the term for interaction in 
the TOIR. We consider a two cases it are ID configuration 
and quasi- ID trap. Let us notice the limiting cases of exis- 
tence of the solution. One of the limiting cases in the region 
of parameters when the scattering length (SL) a is of the or- 
der 10~*cm, and the corresponding equilibrium concentra- 
tion is W^^cm^^. That could be actual in connection with 



1 2711 . Another limiting case is the region of parameters with 
the SL of the order lO^^cm. This case corresponds to con- 
centrations about of 10^^ — lO^^cm"^, which are usually 
dealt with in experiments with BEC. Thus, in order to form 
the conditions for soliton occurrence, the Feshbach reso- 



the development of cooling methods for dense gases 02611 . 



nance (FR) phenomenon should be used 112811 . 112911 . They 
attains the wide-limit SL change in FR experiments, par- 
ticularly the values 10"^ — 10^ Cq (ao-Bohr radius) can be 
reached for magnetically trapped ^^Rb [30]. 

We use in this paper short-range interaction potential 
quantum hydrodynamic equation derived for the system of 
ultracold neutral particles. In connection with this, our at- 
tention should be paid to the fact that an increase in the 
SL can be caused both by a decrease and an increase in 
the depth or width of the interaction potential. Assuming 
that an increase in the SL is caused by a decrease in the 
interaction potential depth, the conditions of existence of 
equations could be considered as fulfilled. 

In a general case, the fact that under the FR condition 
larger values of SL a are attained, can point to the fact 
that a more successive account for the interaction should 
be necessary. 

The processes and effects in the TOIR, along with the 
effects in the spinor BEC lIsTll . magnetically [32], [33], 
[34] and electrical ifssll . [l36ll polarized BEC, can play an 
important role when investigating BEC and interatomic in- 
teraction. 

Our paper is organized as follows. In Sect. 2 we present 
basic equation and describe using model. In Sect. 3 we 
consider a solitons in ID BEC and describe a method of 
getting of solution. We show that with solution is a new 
solution and receive a condition of existence of this solu- 
tion. In Sect. 4 we obtain system of QHD equations for 
the quasi-one dimensional case. In Sect. 5 we investigate 
soliton solution obtained in sect. 3 for quasi-one dimen- 
sional case. In Sect. 6 brief summary of obtained results is 
presented. 



II. MODEL 

To investigate solitons in BEC, we use the set of QHD 
equations up to the TOIR approximation [7]. The calcu- 
lation of the first member in a quantum stress tensor that 
corresponds to the GP equation is fulfilled in [7] under the 
condition that the particles do not interact. A more com- 
plete investigation into the conditions of derivation of the 
GP equation from the MPSE shows that the GP equation 
appears in the first order by the interaction radius (FOIR), 
if the particles are in an arbitrary state that can be simulated 
by a single -particle wave function. Such a state can partic- 
ularly appears as a result of strong interaction between the 
particles that takes place in quantum fluids. 

The QHD equations set for the atoms with a two-particle 
interaction with the potential U{r) and located in external 
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field Vext (fj t) in the TOIR approximation has the form ^ 

dtn{r, t) + 9"(n(r, t)v"{r, t)) = (3) 

and 

mn{r, t)dtV°'(r, t) + ^mn(r, t)5"u^(r, t) 



5"An(r,f) + 

4m 



4m V n(r, t) 



1 



-Tn(r,t)a„n(r,t) - — Ta^aAn (r, t 
Id 



-n(r,f)a"K.t(r,t), 



where 



and 



T 



An 

T 
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dr{r 



dr{r 



dr 



.dU{r) 



dr 



(4) 



(5) 



(6) 



We also have T = —g. Equations Q and (O determine 
the dynamic of concentration of particles n(r, t) and veloc- 
ity field v"(r, t). From equation (|4| we see that dynamics 
of BEC depends on different moments of interaction po- 
tential T, T2. The system of equations (O and is differ 
from (m and (|2]l by existence of one new term. It is a last 
term in left hand side of equation (|4]l. This term appears at 
interaction account up to TOIR approximation. 

In diluted alkali gases, the interaction between parti- 
cles can be considered as scattering. In this case, the 
FOIR interaction constant can be expressed in terms of SL 
T = —4:'jrh^a/m |6]. The second interaction constant T2 
emerges in the TOIR. In a general case, parameter T2 is in- 
dependent of T. In [7], the approximate expression T2 via 
T is considered. We use this expression in our work when 
we investigate the existence region of the soliton solution. 

Considering the dispersion equation for elementary ex- 
citations in BEC accurate to TOIR 



4m" 



+ 



8m 



Tno^2 



m 



which obtained in f?"] we can see that coefficient at fc^ 
could be negative. It is realized at condition T2 < 
—2ff/mno. Consequently, we can expect that value T2 = 
— 2fi,^/mno could play important role at investigation of 
nonlinear processes. 



ni. BRIGHT-LIKE SOLITON IN ID BEC 

In this section, we consider the solitons in the ID BEC. 
For this purpose, we use the perturbative method [i23ll . 
||24f|. Here we present some detail of calculations and de- 
scribe the perturbative method. 

We investigate the case when the stretched variables in- 
clude the expansion parameter in follows combination: 

( = e^/^{z-ut), T = e^/^ut (7) 

where u is the phase velocity of the wave, e- is a small 
nondimension parameter 

An operational relations are arisen from (O 



(8) 



The decomposition of the concentration and velocity 
field involves a small parameter e in the following form: 



n = ng + erii + e 71-2 + ... 



(9) 



V = evi + £^V2 + ... (10) 

Presented in Q equilibrium concentration no is a constant. 
We put expansions (l8]|-([Toll in equations ^ and dill. Then, 
the system of equation is divided into systems of equations 
in different orders on e. 

Equations emerging in the first order by e from the sys- 
tem of equations (|28] | have form 



(11) 



-ud^rii + riod^Vi = 0, 



and lead to the following expression for the phase velocity 



Tno 



m 



(12) 



Square of phase velocity must be positive. Conse- 
quently T is negative, i.e. 



T < 0. 



(13) 



It corresponds to the repulsive SRI. 

Also, from ([TTT l we obtain relation between rii and Vi 
and their derivatives 

OeTli = OcVi. 

U 

Integrating this equation and using a boundary conditions 

ni, f 1 — )• at X — )• ±00 (14) 

we have 

nr. 

(15) 



no 

ni = — vi. 
u 



In the second order by e, from equations (O and dU, we 
derive 

— ud^n2 + udrUi + d^{nQV2 + riiVi) = (16) 



4 



and 



-mu{nod^V2 + riid^vi) + muriQdrVi + mnoVid^Vi 



Am 



Tno(9^n2 + Tuid^Ui + -T2noC?|ni. 

(17) 

In ( [T6l l we can express 712 via f 2 and ni, and put it in 
equation ( fTTl ). Using (fT2l) . we exclude ^2 from the obtained 
equation ( [TTl l. Thus, we obtain an equation which contain 
ni and Vi, only. Using (fTSl l. expressing Vi via ni we get a 
Korteweg-de Vries equation for rii 

drUi + piDTiid^ni + qiod^ni = 0. (18) 

In this equation the coefficients and quy arise in the 
form 

and 



qiD 



2m 



^719X2 



2noT 



(20) 



From equation (fTSl l we can find the solution in the form 
of a solitary wave using transformation r] = ^ — Vt 
and taking into account boundary condition ni = and 

at r/ — )• ±00, we get 



3V 



1 



PlD ,2 
cosn 




(21) 



where V is the velocity of solition propagation to the 
right. From expression pio = 3/2no and solution (1211 
we can find that a perturbation of concentration is posi- 
tive. Consequently, obtained solution is the bright like soli- 
ton (BLS). A width of the soliton is given with formula 
d = 2^yquj/V. BLS exists in the case Quj is positive. 
From condition qiu > ([T3l l we have 



(22) 



Relation (l22l l is fulfil only in the case when T2 is nega- 
tive. In the absence of the second interaction constant T2 
(i. e. in the Gross-Pitaevskii approximation) the relation 
(I22I 1 does not fulfil and, consequently, BLS does not exist. 
From (|22] | we receive that the second interaction constant 
T2 must be negative and it's module must be more than 
Ah'^/mnQ 



IT2I > 



mnQ 



(23) 



Using representation T2 via the s-wave SL a fT*] we get 

T2 = ea^T = -AMa^/m, (24) 

where 9- is a constant, which is determined by an explicit 
form of the interaction potential, > 0, ~ 1 ||7|]. 



From (|22] | and (|24] ) we obtain 

vra^no > 0. 



(25) 



It is the condition of BLS existence. 

Due to used method perturbation rii must be smaller 
than equilibrium concentration tt-q: no >> rii. Here 
we consider the rate ni/no at the centre of soliton at 

cosh( ^17^/27?) = 1: 

ni (centre) 



nn 



Pino 



2V. 



Correspondingly, dimensionless velocity V must be much 
smaller than one. 

At equality in formula (|25] | the BLS has infinite width, 
with the increasing of interaction solitons width becomes 
finite. Formula ( |25] l shows the bound condition for exis- 
tence of soliton. 

Below we consider the same problem for cigar-shaped 
trap. 



IV. THE QUANTUM HYDRODYNAMICS EQUATION IN 
THE CIGAR-SHAPED TRAPS 

Let us to consider the variation of the form of equations 
of QHD in the case of cigar-shaped magnetic traps: 



where loq and Aa;o are angular frequencies in radial and ax- 
ial directions and A is the anisotropy parameter In a quasi- 
ID geometry, the anisotropy parameter in the axially-free 
motion approximation becomes zero A = 0. Thus, the so- 
lution for the radial wave -function appears in the form 



mojQ 



■exp 



(26) 



In the TOIR during the quasi- ID motion of bosons in 
magnetic traps, the GP equation preserves the form but the 
interaction constant [37] changes. A complete 3D particle 
concentration n^(p,z,t) can be presented as product of 
one-dimensional time dependent concentration n{z, t) and 
static radial two-dimensional concentration n{p) 

n^{p, z, t) = n{z, t)n{p), (21) 

where the value n(p) is presented by the formula (l26l l. 

Applying the procedure described in [|37ll . from the set 
of equations (|4| and using the corresponding NLSE, we 
can acquire the system of QHD equations for cigar-shaped 
trap. Here, we describe basic steps of this procedure. Start- 
ing from equation (|4| we get a equation of evolution of a 
following function, which sometimes called wave function 
in medium or order parameter, 



^>(r, t) = n(r, t) ex.p{imd{r, t) /K), 
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where 6 is the potential of velocity field, i.e. v = V9. 
Equation for <l>(r, t) is the NLSE corresponding to system 
of equations (|4]l. Approximately we can present <I>(r, t) in 
the form $(r,t) = ^{p,z,t) = ^{z,t)^{p), where ^{p) 
is the wave function of the ground state of harmonic oscil- 
lator and the square of module of ^{p) presented by for- 
mula (l26l l. Since, we get a NLSE for <I>(2;, t). This equation 
describes the evolution of BEC in quasi-one dimensional 
trap. From obtained NLSE we derive the system of QHD 
equations for quasi- ID trap. In the results we have 



and 



dtn{z,t) + dz{n{z,t)v{z,t)) =0 



mn{z, t)dtv{z, t) H — mn{z, t)dzv'^{z, t) 



4m 4m n{z,t) 

+ain{z, t)dzn{z, t) + a2n{z, t)dln{z, t) 

-a2{dMz,t))d,n^{z,t)-a2 ^ " / \ '' = 0. 



(28) 



The following parameters appear in equation 



evi + e V2 + 



(31) 



Presented in (|30] | equilibrium concentration Uq is a con- 
stant. Equations emerging in the first order by e from the 
set of equations (|28] | lead to the following expression for 
the phase velocity u: 



u 



m 



m\ ZTrn, 2 \ ttH J J 



It is evident from ( |32] | that the wave can exist under the 
condition ai > 0. The obtained condition means that in 
the repulsive forces should act in the case under consider- 
ation under the condition that the contribution of terms in 
FOIR prevails over the TOIR terms. 

Relationship (l32l l is the analog of the dispersion depen- 
dence. The use of scaling ( [29l l leads to simplifying the dis- 
persion relationship compared with the case when we con- 
sider small perturbations proportional to exp(—iujt + ikz). 
In the latter case, we obtained the dispersion relation uj{k) 
form the set of equations (|28] | in the form: 



U! 



4m 



16m ttH 



and 



1 mwo Stt / muo 

«! = — T 19 

2 irh 2 "I TTfi 



3 „ muo 
«2 = -—^2- 



16 nh 

The form of nonlinear terms that describe the interaction 
in equation (l28l l differs from corresponding terms in (@). 
This leads to varying the form of solutions and conditions 
of their existence in a quasi- ID geometry compared with a 
ID case. 



V. THE SMALL AMPLITUDE SOLITONS IN QUASI-ID 
BEC 

In this section, we consider the solitons in the BEC 
for the case of small nonlinearity taking into account the 
TOIR. For this purpose, we use the perturbative method 

We investigate the case when the stretched variables in- 
clude the expansion parameter in follows combination: 



{z — Ut), T 



(29) 



where u is the phase velocity of the wave. 

The decomposition of the concentration and velocity 
field involves a small parameter e in the following form: 



n 



no + erii + e n2 + 



(30) 



Thus, relationship (|32] | corresponds to the phonon part of 
the dispersion dependence (|33] |. 

From the second-order set of equations (|2Ft by e we find 
that the concentration rii satisfies the Korteweg-de Vries 
equation: 

d^rix +pnid^ni + qd^Ui = 0, (34) 

where 

3 

P=l^, (35) 



2no' 

4mnoa2 — 
SmnoOi 



(36) 



When we derive this equation, we used boundary condi- 
tions 111 = and Vi = at ^ — )■ ±00. Using transforma- 
tion Tj = ^ — Vt, and taking into account boundary condi- 
tion 77-1 = and d^rii = at — t- ±00, we can obtain the 
solution in the form of a solitary wave from equation (|34) 



ni 



31/ 



1 



^ cosh^l 



(37) 



where V is the velocity of solition propagation to the right. 
Sign of perturbation is determined by the sign of p. From 
formulas (|35] |. (|36] | we can see the quantity p is positive. 
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Consequently, obtained solution is the soliton of compres- 
sion or bright like soliton solution, by analogy with well- 
known bright soliton in BEC Hill, 111, S. As it will 
be shown below, this solution exists only when taking into 
account the TOIR. 

Let us pass on to a detail consideration of the conditions 
of existence of the solution ( l37l) . The solution (l37l l of the 
equation (l34l l exists as the conditions q > and ai > are 
fulfilled. We start with consideration the condition q > 0. 
As Oi > 0, then to fulfill the condition g > we need 
— + 4mnoa2 > 0. 

In the case when 02 is vanish (i.e. in FOIR approxima- 
tion) the solution (l37l) is not exist. It means, that solution 
arises in TOIR approximation which developed in [7]. One 
of the condition of existence of the solution (l37l l is: 

a2 > . (38) 

Consequently, for the second interaction constant T2 we 
obtain: 

T2 < -— . (39) 

Using relation (|24] | we can make estimation for corre- 
sponding SL. It is useful to present the value of possi- 
ble SL in the terms of space parameter of the trap a± = 
^hjrabjf^. From conditions 
SL a appear: 



the conditions for 



a > 



39no 



In addition, from > and (l24l l. we obtain 

a 

a < 



56 



(40) 



(41) 



Using equation (l24ll . the particle concentration n can be 
presented in the form 



n 



where 



no + 2Vnoe • sech^ ( W~7^ ) ' 



39mnQUjQa^h — 



(42) 



(43) 



The soliton width d arises in the form d = 2^J q^ jV . 
The numerical analysis of formula (|43^ is presented in Fig. 
[TJIJ] It is evident from Fig. [HE] that there is a narrow in- 
terval of the SL values, for which the solution (|42] |. (|43^ 
exists. The dependence of the soliton width on the SL 
is resonant-shaped. The resonant value of the SL de- 
pends on the particle concentration no and trap parameter 
The values of become lower at increasing of equilibrium 
concentration no. The SL reaches 0.1 nm at the concen- 



As the concentration de- 

-3 



tration of the order 10 cmT 

creases to values 10^^-10^"* cm^'^, which are usually used 
in BEC experiments, the SL increases to the values of the 
order of 1 fim. Such values of the SL can be attained when 
using the FR. 




-6.0 -5.5 -5.0 -4.5 
In(fl), [a\ = cm 




-5.354 -5.352 -5.350 -5.348 -5.346 -5.344 
ln(a), [«] = cm 

FIG. 1. The dependence of the soliton width d on the scatter- 
ing length a at fixed parameter of the trap — ^Jh/mujf) = 
10~^cm and equilibrium concentration no = Kf'cm^^ and at 
V — 1, 9 — 1. On Fig. la we can see what width of soliton is the 
positive in small range of the values of the SL. On Fig. lb the area 
of the resonance is presented more detailed than on Fig. la. 



6.0 




[a^] = cw -5,62 



-5.60 



FIG. 2. (Color online) The dependence of soliton width d on ra- 
dial parameter of the trap aj^ and the nonperturbative concentra- 
tion of the particles no at the fixed scattering length a — lO^^cm, 
and atV = I, 6 = 1. The soliton width d is positive in two area. 
But in area with smaller value of the particles concentration no 
the square of phase velocity is negative. Thus, the solution exist 
in area of bigger concentrations. 
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VI. CONCLUSION 

In this article, we showed that at a more exact account- 
ing of the interaction, specifically, taking into account the 
TOIR, a new type of solitons emerges in the BEC. We 
also studied the conditions of existence of such a solution. 
For this problem solving we used the set of QHD equa- 
tions where the interactions included up to TOIR approx- 
imation. The TOIR approximation is an example of the 
nonlocal interaction. The GP approximation gives us the 
force density in the right hand side of the Euler equation 
F = —gVn?/2 = TVn^/2. It is corresponds to the first 
order on interaction radius. The interaction including up to 
TOIR approximation gives us the second tern in the force 
field F = TVnV2 + T2VAnVl6. The new term contain 
the third spatial derivative of the concentration square and 
the new interaction constant. For obtained results analysis 
and estimation we used approximate estimation of T2 and 
its approximate connection with the T or SL a. 

We found that BLS (soliton of compression) exists in 
ID case (one dimensional propagation in three dimensional 
medium) in the case of strong enough repulsive interaction. 
BLS appearance strongly connects with the second interac- 
tion constant T2. If we consider the interaction in the first 
order by the interaction radius there is no BLS. We also 
studied the BLS behavior in the case of quasi- ID trap and 
describe contribution of external fields on BLS amplitude 
and width. 

In a general case, the second interaction constant T2, 
which appears in the TOIR, is independent of T and, con- 
sequently, of SL a. Thus, the relationships obtained in this 
article (|35] |. (|36] |. (|37] | can be used for an independent ex- 
perimental determination of T2. In this case, parameter 
h^/mn^o can be used for the qualitative evaluation of the 
second interaction constant T2. 

Thus, in this paper we showed that new physical effects 
appear at account of interaction up to TOIR approxima- 
tion. The processes and effects in the TOIR approximation, 
along with the effects in the spinor and polarized BEC, can 
play an important role at investigation of BEC and inter- 
atomic interaction. 
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